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A Kubo-Greenwood-like equation for the Gilbert damping parameter a is presented that is based 
on the linear response formalism. Its implementation using the fully relativistic Korringa-Kohn- 
Rostoker (KKR) band structure method in combination with Coherent Potential Approximation 
(CPA) alloy theory allows it to be applied to a wide range of situations. This is demonstrated with 
results obtained for the bcc alloy system Fe^Coi^ as well as for a series of alloys of permalloy with 
5d transition metals. To account for the thermal displacements of atoms as a scattering mechanism, 
an alloy-analogy model is introduced. The corresponding calculations for Ni correctly describe the 
rapid change of a when small amounts of substitutional Cu are introduced. 

PACS numbers: Valid PACS appear here 



I. INTRODUCTION 

The magnetization dynamics that is relevant for the 
performance of any type of magnetic device is in gen- 
eral governed by damping. In most cases the magneti- 
zation dynamics can be modeled successfully by means 
of the Landau-Lifshitz-Gilbert (LLG) equation [l[ that 
accounts for damping in a phenomenological way. The 
possibility to calculate the corresponding damping pa- 
rameter from first principles would open the perspective 
of optimizing materials for devices and has therefore mo- 
tivated extensive theoretical work in the past. This led 
among others to Kambersky's breathing Fermi surface 
(BFS) [2] and torque-correlation model (TCM) 0, that 
in principle provide a firm basis for numerical investi- 
gations based on electronic structure calculations [J, Q . 
The spin-orbit coupling that is seen as a key factor in 
transferring energy from the magnetization to the elec- 
tronic degrees of freedom is explicitly included in these 
models. Most ab-initio results have been obtained for the 
BFS model though the torque-correlation model makes 
fewer approximations [HQ. In particular, it in principle 
describes the physical processes responsible for Gilbert 
damping over a wide range of temperatures as well as 
chemical (alloy) disorder. However, in practice, like many 
other models it depends on a relaxation time parame- 
ter r that describes the rate of transfer due to the vari- 
ous types of possible scattering mechanisms. This weak 
point could be removed recently by Brataas et al. Q 
who described the Gilbert damping by means of scatter- 
ing theory. This development supplied the formal basis 
for the first parameter-free investigations on disordered 
alloys for which the dominant scattering mechanism is 
potential scattering caused by chemical disorder Q. 

As pointed out by Brataas et al. Q , their approach is 
completely equivalent to a formulation in terms of the lin- 
ear response or Kubo formalism. The latter route is taken 



in this communication that presents a Kubo-Greenwood- 
like expression for the Gilbert damping parameter. Ap- 
plication of the scheme to disordered alloys demonstrates 
that this approach is indeed fully equivalent to the scat- 
tering theory formulation of Brataas et al. @- In addi- 
tion a scheme is introduced to deal with the temperature 
dependence of the Gilbert damping parameter. 

Following Brataas et al. 0, the starting point of our 
scheme is the Landau-Lifshitz-Gilbert (LLG) equation 
for the time derivative of the magnetization M: 
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where M s is the saturation magnetization, 7 the gyro- 
magnetic ratio and G the Gilbert damping tensor. Ac- 
cordingly, the time derivative of the magnetic energy is 
given by: 
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in terms of the normalized magnetization m = M/M s . 
On the other hand the energy dissipation of the electronic 

system E^is = (^Sr^j ^ s determined by the underlying 

Hamiltonian H(t). Expanding the normalized magne- 
tization to(t), that determines the time dependence of 
H(t) about its equilibrium value, m(r) = too + u(t), one 
has: 



H = H (m ) + 



d_ 

1 du n 



(3) 



Using the linear response formalism, E^is can be written 
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where .Bp is the Fermi energy and the sums run over all 
eigenstates a of the system. Identifying -E mag = E&s, 
one gets an explicit expression for the Gilbert damping 
tensor G or equivalently for the damping parameter a = 
G/hMs): 
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(5) 



An efficient way to deal with Eq. ([5]) is achieved by ex- 
pressing the sum over the eigenstates by means of the 
retarded single-particle Green's function lmG + (E F ) — 
-7r^ a \if> a ){ij; a \S{E F — E a ). This leads for the parame- 
ter a to a Kubo-Greenwood-like equation: 



--^Trace ( |^-Im G+(£ F )f^Im G+(E F ) 
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with (...) c indicating a configurational average in case of 
a disordered system (see below). Identifying dH/du^ 
with the magnetic torque T M this expression obviously 
gives the parameter a in terms of a torque-torque corre- 
lation function. However, in contrast to the conventional 
TCM the electronic structure is not represented in terms 
of Bloch states but using the retarded electronic Green 
function giving the present approach much more flexibil- 
ity. As it corresponds one-to-one to the standard Kubo- 
Greenwood equation for the electrical conductivity, the 
techniques developed to calculate conductivities can be 
straightforwardly adopted to evaluate Eq. ([5]). 

The most reliable way to account for spin-orbit cou- 
pling as the source of Gilbert damping is to evaluate 
Eq. (J6j) using a fully relativistic Hamiltonian within the 
framework of local spin density formalism (LSDA) Q : 



H = cap + (3mc 2 + V(r) + (3amB(r) . 



(7) 



Here on and j3 are the standard Dirac matrices and p is 
the relativistic momentum operator (To| . The functions 
V and B are the spin-averaged and spin-dependent parts 
respectively of the LSDA potential. Eq. ([7]) implies for 
the magnetic torque Tp occurring in Eq. © the expres- 
sion: 



Tfi = -£-H = /3Ba n . 



(8) 



The Green's function G + in Eq. ([5]) can be obtained in a 
very efficient way by using the spin-polarized relativistic 



version of multiple scattering theory [TJ] that allows us to 
treat magnetic solids: 



G+(r n ,r m ',E) = £ ZJ£(F„, E) (E) Z™ x (f m E) 

AA' 

- Y, z l <f< . E ) j a< X (f> > E) S nm .(9) 

A 

Here coordinates r n referring to the center of cell n 
have been used with |r<| = min(\r n \, \r n '\) and |r>| = 
maa;(|r n |, \f n '\). The four component wave functions 
Z^(f,E) (Jj^(f,E)) are regular (irregular) solutions to 
the single-site Dirac equation for site n and Tp^,{E) is 
the so-called scattering path operator that transfers an 
electronic wave coming in at site m into a wave going 
out from site n with all possible intermediate scattering 
events accounted for coherently. 

Using matrix notation, this leads to the following ex- 
pression for the damping parameter: 
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with the g- factor 2(1 + /i or b/ 'fi spin) m terms of the spin 
and orbital moments, ii sp i n and fiorb, respectively, the 
total magnetic moment fi tot = \i sp i n + n or b, and Tjfy, — 
??( T AA' ~ t a?a) anc ^ * ne ener gy argument E F omitted. 
The matrix elements of the torque operator T ra/i are iden- 
tical to those occurring in the context of exchange cou- 
pling fllj | and can be expressed in terms of the spin- 
dependent part B of the electronic potential with matrix 
elements: 



rpUfJ, 

1 A'A 



d\ (r) \j3a^B xc (f)] Zftf) . (11) 



As indicated above, the expressions in Eqs. ([6]) - (fTTj) 
can be applied straightforwardly to disordered alloys. 
In this case the brackets (...) c indicate the necessary 
configurational average. This can be done by describ- 
ing in a first step the underlying electronic structure 
(for T — K) on the basis of the Coherent Potential 
Approximation (CPA) alloy theory. In the next step 
the configurational average in Eq. ([6]) is taken follow- 
ing the scheme worked out by Butler [12j when dealing 
with the electrical conducting at T = K or residual 
resistivity respectively, of disordered alloys. This im- 
plies in particular that so-called vertex corrections of the 
type (T M ImG+T„ImG+) c - (T u lmG+) c (TJmG+) c that 
account for scattering-in processes in the language of the 
Boltzmann transport formalism are properly accounted 
for. 

Thermal vibrations as a source of electron scattering 
can in principle be accounted for by a generalization of 
Eqs. - (fTTj) to finite temperatures and by includ- 
ing the electron-phonon self-energy S e /_p^ when calcu- 
lating the Greens function G + . Here we restrict our- 
selves to elastic scattering processes by using a quasi- 
static representation of the thermal displacements of the 
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FIG. 1: Gilbert damping parameter for bcc Fe x Cor-a; as a 
function of Co concentration: full circles - the present results 
within CPA, empty circles - within non-local CPA (NL CPA), 
and full diamonds - experimental data by Oogane [T3 ]. 

atoms from their equilibrium positions. We introduce 
an alloy-analogy model to average over a discrete set 
of displacements that is chosen to reproduce the ther- 
mal root mean square average displacement y/ (u 2 )t for 
a given temperature T. This was chosen according to 

<Ot = l^m^i^ff 1 + W -th *(0 D /T) the De- 
bye function, h the Planck constant, k the Boltzmann 
constant and the Debye temperature [13j. Ignoring 
the zero temperature term 1/4 and assuming a frozen 
potential for the atoms, the situation can be dealt with 
in full analogy to the treatment of disordered alloys de- 
scribed above. 

The approach described above has been applied to 
the ferromagnetic 3d-transition metal alloy systems bcc 
Fe^Coi-j;, fee Fea.Nii._2. and fee Co a Nii_3;. Fig. [TJshows 
as an example results for bcc Fe^Coi-a; for x < 0.7. The 
calculated damping parameter a(x) for T = K is found 
in very good agreement with the results based on the 
scattering theory approach [8] demonstrating numerically 
the equivalence of the two approaches. An indispensable 
requirement to achieve this agreement is to include the 
vertex corrections mentioned above. In fact, ignoring 
them leads in some cases to completely unphysical re- 
sults. To check the reliability of the standard CPA, that 
implies a single-site approximation when performing the 
configurational average, we performed calculations on the 
basis of the non-local CPA 15]. In this case four atom 
cluster have been used leading - apart from the very di- 
lute case - practically to the same results as the CPA. As 
found before for fee Fe^Mi-^ || the theoretical results 
for a reproduce the concentration dependence of the ex- 
perimental data quite well but are found too low (see 
below). As suggested by Eq. (ITU)) the variation of a(x) 
with concentration x may reflect to some extent the vari- 
ation of the average magnetic moment fj, tot of alloy. As 
the moments as well as the spin-orbit coupling strength 
of Fe and Co don't differ too much, the variation of a(x) 
should be determined in the concentrated regime primar- 
ily by the electronic structure at the Fermi energy Ep- 



As Fig. [TJshows, there is indeed a close correlation of the 
density of states n(Ep) that may be seen as a measure 
for the available relaxation channels. 

While the scattering and linear response approach are 
completely equivalent when dealing with bulk alloys the 
latter allows us to perform the necessary configuration 
averaging in a much more efficient way. This allows us to 
study with moderate effort the influence of varying the 
alloy composition on the damping parameter a. Corre- 
sponding work has been done in particular using permal- 
loy as a starting material and adding transition metals 
(TM) [l6| or rare earth metals [l?} • Fig. [3] (top) shows re- 
sults obtained by substituting Fe and Ni atoms in permal- 
loy by 5d TMs. As found by experiment [l6| a increases 
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FIG. 2: Top: Change of the Gilbert damping parameter Aa 
w.r.t. permalloy (Py) for various Py/5d TM systems as a func- 
tion of 5d TM concentration; Middle: Gilbert damping pa- 
rameter a for Py/5d TM systems with 10 % 5d TM content 
in comparison with experiment [Ji|; Bottom: spin magnetic 
moment m 5 a p in and density of states u(Ef) at the Fermi en- 
ergy of the 5d component in Py/5d TM systems with 10 % 
5d TM content. 

in all cases nearly linearly with the 5d TM content. The 
total damping for 10 % 5d TM content shown in the 
middle panel of Fig. [5] varies roughly parabolically over 
the 5d TM series. In contrast to the Fe^Co!-^ alloys 
considered above, there is now an S-like variation of the 
moments over the series (Fig. [5J bottom), char- 

acteristic of 5d impurities in the pure hosts Fe and Ni 
[Isl [l9| . In spite of this behaviour of fJ.^p in the variation 
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FIG. 3: Temperature variation of Gilbert damping of pure 
Ni and Ni with Cu impurities: present theoretical results vs 
experiment [20] 

of a(x) seems again to be correlated with the density of 
states n 5d (E F ) (Fig. [5] bottom). Again the trend of the 
experimental data is well reproduced by the theoretical 
ones that are however somewhat too low. 

One of the possible reasons for the discrepancy of the 
theoretical and experimental results shown in Figs. [1] and 
[2] might be the neglect of the influence of finite temper- 
atures. This can be incorporated as indicated above by 
accounting for the thermal displacement of the atoms in 
a quasi-static way and performing a configurational av- 
erage over the displacements using the CPA. This leads 
even for pure systems to a scattering mechanism and this 
way to a finite value for a. Corresponding results for pure 
Ni are given in Fig. [3] that show in full accordance with 
experiment a rapid decrease of a with increasing tem- 
perature until a regime with a weak variation of a with 
T is reached. This behavior is commonly interpreted as 
a transition from conductivity-like to resistivity-like be- 
haviour reflecting the dominance of intra- and inter-band 
transition, respectively Q , that is related to the increase 
of the broadening of electron energy bands caused by the 
increase of scattering events with temperature. Adding 
only less than 1 at. % Cu to Ni, the conductivity-like 
behavior at low temperatures is strongly reduced while 



the high temperature behavior is hardly changed. A fur- 
ther increase of the Cu content leads to the impurity- 
scattering processes responsible for the band broaden- 
ing dominating a. This effect completely suppresses the 
conductivity-like behavior in the low-temperature regime 
because of the increase of scattering events due to chem- 
ical disorder. Again this is fully in line with the experi- 
mental data, providing a straightforward explanation for 
their peculiar variation with temperature and composi- 
tion. 

From the results obtained for Ni one may conclude that 
thermal lattice displacements are only partly responsible 
for the finding that the damping parameters obtained 
for Py doped with the 5d TM series, and Fe x Coi_2; are 
somewhat low compared with experiment. This indicates 
that additional relaxation mechanisms like magnon scat- 
tering contribute. Again, these can be included at least 
in a quasi-static way by adopting the point of view of a 
disordered local moment picture. This implies scatter- 
ing due to random temperature-dependent fluctuations 
of the spin moments that can also be dealt with using 
the CPA. 

In summary, a formulation for the Gilbert damping pa- 
rameter a in terms of a torque-torque-correlation func- 
tion was derived that led to a Kubo-Greenwood-like 
equation. The scheme was implemented using the fully 
relativistic KKR band structure method in combination 
with the CPA alloy theory. This allows us to account for 
various types of scattering mechanisms in a parameter- 
free way. Corresponding applications to disordered tran- 
sition metal alloys led to very good agreement with re- 
sults based on the scattering theory approach of Brataas 
et al. demonstrating the equivalence of both approaches. 
The flexibility and numerical efficiency of the present 
scheme was demonstrated by a study on a series of 
permalloy-5d TM systems. To investigate the influence of 
finite temperatures on a, a so-called alloy-analogy model 
was introduced that deals with the thermal displacement 
of atoms in a quasi-static manner. Applications to pure 
Ni gave results in very good agreement with experiment 
and in particular reproduced the dramatic change of a 
when Cu is added to Ni. 
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